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Abstract 

We derive the forward photoproduction amphtude for the diffractive jp l~^l~p reaction in 
the momentum space, within the formaUsm of factorization. Predictions for the "fp — >■ l~^l~p 
reaction are given using unintegrated gluon distribution from the Uterature. We calculate the 
total cross section as a function of photon-proton center of mass energy and the invariant mass 
distribution of the lepton pair. We also discuss whether the production of timelike virtual photons 
can be approximated by continuing to the spacelike domain < 0. The present calculation 
provides an input for future predictions for exclusive hadroproduction in the pp — )• p l^l~p reaction. 
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I. INTRODUCTION 



Measuring absolutely normalized cross sections at the LHC is of great importance for the 
high-energy physics community. This requires having a well understood luminosity monitor. 
Following the pioneering work the QED process pp — )■ pl^l~p via photon-photon fusion is 
often discussed as a process which can be used for measuring the luminosity at the LHC 0- 
3]. It is therefore very important to estimate other non-QED contributions to exclusive 
production. One possible source of dileptons is the exclusive production of vector mesons or 
Z-bosons (see e.g. [5-7]). The dilepton pairs originating from these processes however have 
invariant masses close to the mass of the decaying state. In Figure [1] we show an exclusive 
diffractive mechanism which produces continuum dilepton pairs, and hence may compete 
with the standard QED process. In this reaction the coupling of the photon to the proton 
is known and can be expressed in terms of the nucleon electromagnetic form factors. At the 
small transferred momenta (ti or ^2) of relevance, it is sufficient, in the high energy limit, 
to include the Dirac electromagnetic form factor. Besides its role as a possible background 
to electromagnetic lepton pair production, the 7p — ?• 7*p amplitude may contain interesting 
information on the small-x gluon distribution in the nucleon. 




FIG. 1: An example of the non-QED mechanism for the prodution of opposite charge leptons in 
the pp — 7> ppl^l~ reaction. 

In the present work we shall concentrate on the photon-pomeron subprocess. In Fig|2]we 
show a QCD mechanism, where the photon splits into a quark-antiquark pair which interacts 
with the proton through the exchange of an off-diagonal QCD gluon ladder. In principle this 
process could have been studied at HERA. In Figj2] the incoming photon is spacelike, (or 
quasireal) but the outgoing photon is timelike, i.e. its virtuality > 0. This process is often 
called timelike Compton scattering (TCS) in the literature, although the specific mechanism 
considered by us is maybe better termed a QCD version of (virtual) Delbriick scattering. A 
coUinear factorisation treatment of timelike Compton scattering in terms of the nucleon's 
skewed (mainly quark-) distributions can be found in [9|. This approach is most relevant 
for lower center-of-mass energies. We will restrict ourselves to high energies, where the 
t-channel exchange is dominated by gluons, and choose a /c^-factorization formalism very 
similar to the one used in diffractive vector meson production [13, U . 



The TCS cross section has also been evaluated in a color-dipole model with a saturation- 
idea inspired dipole-nucleon cross section [l^. However there both incoming and outgoing 
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photons were assumed to be spacelike. A recent estimate of high energy cross-sections in 
leading order collinear factorization, without explicit gluons, is found in [isl]. 

In this work we present the momentum-space formulation of timelike Compton scattering 
at small-x, taking due account of the timelike nature of the final state photon. 

This paper is organized as follows. In the next section, we present the formalism which is 
used in our calculations. In Sec. Ill we present the main results for the 7p — i- l^l~p reaction. 
Finally, in the last section we summarize our results and show further perspectives. 



II. FORMALISM 




FIG. 2: The diagram for the production of virtual timelike photons. 



The photoproduction amplitude will be the major building block for our prediction of 
exclusive dilepton pair production. The amplitude for the reaction is shown schematically 
in Figj2l In the diagram, we distinguish three stages of the process: first the incoming real 
photon fluctuates into a quark-antiquark pair, then a gluon ladder is exchanged between the 
qq pair and the proton and finally the qq pair recombines to form a virtual photon which 
subsequently decays into a lepton-antilepton pair. 

The amplitude of the subprocess jp — )■ 'j*{q'^)p is a sum of the contributions for a given 
flavour / of quarks in the loop. 

M{-fp -> 7*(g')p) = J2^f^^P ^ (2.1) 

/ 

Here the initial state photon is real, and hence transversely polarized. We take only the 
dominant s-channel helicity conserving contribution into account, and suppress helicities 
of photons/protons in our notation. The calculation of the amplitude follows the same 
procedure as for the exclusive production of vector mesons, which is explained in great detail 
in Ivanov's thesis [lo|. The main difference is that the final state light-cone wavefunction 
is replaced by a free quark propagator times the QED-spinor structure for the qq — )■ 7* 
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transition. The forward 7J9 — )■ 7*p amplitude for a given flavour contribution can then be 
written as: 



Jo Jo 



[k\ + m'j — z{l — z)q'^ — ie] 



dz r ..o Afiz,kl,W^) 



z{l-z) 



Tidk 



— 7-i — f — —le 

2(1 — Z) ^ 



, (2.2) 



where the exphcit form of Af{z, k\,W'^) will be discussed below, agm is the QED fine- 
structure constant; e/ = | for u,c,t and e/ = — | for d,s,b is the quark charge. The 
transverse momentum squared of (anti-)quarks is denoted by A;^, their longitudinal momen- 
tum fractions are z and 1 — z, respectively. The running coupling as enters at the scale 
= maxjfi;^, fc^ + m'j}, where rrif is the quark mass for flavor /. Now notice, that the 
invariant mass of the qq pair is given by 



„ k^ + ml 



(2.3) 



z{l - z) ' 

so that the second line of Eq. fl2.2p suggests a change of variables from (z, k\) — )■ (M^, k\)\ 

dzdk\ dM'^ dk\ 

z{i-z) 'W ~T 

Using the z -^r^ 1 — z symmetry we could restrict ourselves to < 2; < 1/2, so that from 



z = z{M\ki) 



1-4- 



/ 



M2 



(2.5) 



we obtain the jacobian factor 



J 



1-4 



/c]_ + m^ 
M2 



which introduces an integrable singularity (see e.g. Eq. 
transformed as 



(2.6) 



The integration domain is 



{0 < 2 < 1/2} X {0 < A;i < 00} — > {Am] < < 00} x {0 < /c^ < {M"^ - Am]) / A] . (2.7) 



Finally, we can cast the amplitude in the form 



-M/(7P^7*(g» = W^'l6vr2«emeJ ■ / 

J4i 



4mj M'^ -q^ - ie 



dM^ 



(2.8) 



Here ajiW'^ , M'^) is related to the diffractive amplitude for the 7p — )■ qqp transition [14 
however with the spinorial contractions from the final state performed. For the lack of a 
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better name we will refer to it as the spectral distribution or spectral density. Its imaginary 
part is given by the integral: 



J 



'^mAf{z{M%ki),ki,W'). 



with 



TX 



+A^f{z,kl)Wij{kl,K' 
where the auxiliary functions ^1,0/5 W^i,o/ can be taken from Ref. [6j 



(2.9) 



(2.10) 



Aof{z,kl) = mj, A^f{z,kl) = [z^ + {l~zf] 



Wof{ki,K') 



k^ + ' 



W^i/(A;i,«:'^ 



/ 



2/t 



1 + 



k'j_~ mj — 



{k1_ — m'j — H?Y + Amj-k^ 



(2.11) 



Given the relation to the diffractive •yp — )■ qqp amplitude we feel justified to obtain the real 
part from the standard derivative form of the dispersion relation: 

neMz, kl, = ^mAflz, kl, W^) ■ tau {f-}2i^^A^l^) . (2,12) 

The function J^{x, h?) is the unintegrated gluon distribution of the proton, which at large 
values of gluon transverse momenta can be expressed in terms of the coUinear gluon distri- 
bution as 



dG{x^ 
d log 



(2.13) 



In the present analysis we use an unintegrated gluon distribution obtained by fitting to 
the structure function data measured at HERA 15|. Following 10|, lll|], we correct for 



skewedness effects by taking the unintegrated at x^s = Cskewed ^ with Cskewed = 0.41. To 
calculate the integral in Eg. (12. 8 p for > 4m^ we use the Plemelj-Sokhocki formula: 



= PV- + m6{x) 
X — ie X 



(2.14) 



where PV denotes the principal value integral. We can finally represent real and imaginary 
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part of the forward TCS amplitude as 



f ^ 
2 nj2) 



^ J Am?, 



f ^ 

2 /\#2\ 



+e{q'-Am}) (PV r ^M^ ^"7j^^'f ^ - vr 53m aKW^^ g^)) | 
V J4m2 M — q /J 



(2.15) 

Two comments are in order on the final form of the TCS amplitude fl2.15p : firstly, the "ivr"- 
terms from the decomposition f l2.14p lead to a quite nontrivial structure of the amplitude 
in terms of the spectral distribution aj. In particular, the TCS amplitude will not be 
a monotonous function of g^. Secondly, we should remember that these "ivr" -terms derive 
from the cut through the qq pair which exists in the perturbative amplitude, when > Am'j. 
Clearly in the real hadronic world, there are no cuts due to quarks going on-shell, and our 
amplitude must, as usual, be interpreted in a parton-hadron duahty sense. 

As a last step we must extend our amplitude to finite momentum transfers. For simplicity 
we will assume it to have the following factorized form 

Mfi^p l*{q^)p-,t) =Mf{ip ^ l*{q^)p) exp[5t] , (2.16) 

which should be sufficiently accurate for small t, within the diffraction cone. The total cross 
section for the 7p — )■ 7*p process can then be obtained as 

a{-fp ^ 7 (g )P) = .a u — ' (2.17) 

where for a first estimation we shall take B = A GeV~^. 

The invariant mass distribution of dileptons for the 7p — )■ l^l~p reaction which can be 
accessed in experiment is given by: 



da 

dq^ ^ Snq 



,{1P ^ in~p) = ■ a{^p -> 7*(g')p), (2.18) 



where means either e'^e or /x^/x . This simple formula applies only when g^ ^ mf. 



III. RESULTS 



Here we present predictions for the jp — )■ l^l~p reaction. In Figj3]we show spectral density 
as a function of (invariant mass of the gg pair) for different 7p center of mass energies: 
W = 100, 500, 1000 GeV and for u, d, s, c quarks separately. The energy dependence of the 
spectral density derives from the a;efr- dependence of the unintegrated gluon density. This 
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FIG. 3: Spectral density (aj) for different flavours: solid line for u,d, dashed line for s, dotted line 
for c at W =100, 500, 1000 GeV. In the left row we show the imaginary part '^maj, while in the 
right row, we show the real part Ikeaj. 

is why we observe a growth of the spectral density with energy. For large invariant mass of 
the qq pair the spectral density tends to zero, which ensures the convergence of the integrals 
in ^IM- 

In Fig. Hlwe show the invariant mass distribution da/dq^ of dileptons in the 7p l^l^p 
reaction as a function of photon-proton center-of-mass energy at fixed values of the invariant 
masses of the dilepton pairs {q"^)- Here all flavours are included in the amplitude. In general, 
the higher invariant mass, the faster growth with the photon-proton energy. This points to 
the fact that the unintegrated glue is probed at on average harder scales, where it has a 
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faster x-dependence. 




FIG. 4: The cross section for ^ l^l p as a function of center-of-mass photon-proton energy for 
fixed values of the dilepton invariant mass. 




FIG. 5: atothp ^ l+l-p] for (^Ln = 0-l> 0.5, 1.0, 1.5 GeV^. 

In Fig. |5]we show cross section integrated over dilepton invariant mass 

poo 
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These cross sections are by a factor of about 5 larger than those in [12| . 

In Figj6] we show the invariant mass distribution of dileptons da / dq^ as a function of 
dilepton invariant mass at fixed values of 7p energy. An interesting feature of the invariant 
mass distribution is a cusp at ~ 9GeV^. Notice that this is the vicinity of q^ ~ 4m^. 
Indeed the cusp is caused by the cc contribution to the amplitude which changes sign in 
this region. This is a unique feature of the structure of the amplitude (12.1 5p with timelike 
final state photons. Here one should however remember the caveat on the absence of quark 
thresholds, most optimistically one may hope that such a structure survives in the vicinity 
of q^ ~ 4mf). Finally, it is interesting to investigate how well the dilepton mass spectrum 




FIG. 6: Distribution in for W = 100, 500, 1000 GeV. 

can be calculated from the amplitude for production of spacelike photons in the final state. 
In this case one would replace Eq. ( 12.151) by the straightforward 

MJ-'^^^^P ^ l*{q')p) = W leTT^aemeJ ■ H d.M\ (3.2) 

J4m2 ^Vl -j- q 

Notice that in distinction to the timelike case, this is a monotonous function of q^. In Fig. [7] 
we show the ratio of da /dq^ for the (correct) timelike photons and the spacelike prescription 
as a function of q^ at various fixed energies. We observe that, as expected, the spacelike 
prescription does not reproduce the structure present in the timelike amplitude. At large q^ 
it gives a fairly reasonable description of the g^-dependence, but is not able to reproduce 
the correct timelike results. The cross section for timelike photons is bigger by a factor of 
3-4 compared to the spacelike photon prescription. 



IV. CONCLUSIONS 



We have derived the amplitude for the exclusive diffractive photoproduction of lepton 
pairs in the /c^-factorisation approach in the momentum space. We have discussed several 
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FIG. 7: The ratio of the cross section -;- as a function of o . 

spacelike 



details of the formalism as well as differences compared to the existing calculation in the 
literature which ignored the fact that the " produced" photons are timehke. 

We have calculated the cross sections as a function of photon-proton center of mass energy 
as well as a function of dilepton invariant mass. We have demonstrated how important is 
the inclusion of correct dynamics (timelike outgoing photons instead of spacelike outgoing 
photons). As a consequence the cross sections obtained here are significantly larger than 
those obtained in the literature. We furthermore observed an interesting structure in the 
invariant mass distribution of dilcptons around ~ 9 GeV^. 

The amplitude for the 7p l^l~p is the main ingredient of the diffractive amplitude 
for the pp — )■ l^l^pp process. In the case of hadroproduction of dileptons the diffractive 
mechanism constitutes a background to the purely electromagnetic (photon-photon fusion) 
production of dileptons. The latter process was sugested in the literature as a luminosity 
monitor the LHC studies. How important is this background for differential distributions in 
the four-body pp l^l^pp reaction will be studied in detail elsewhere. 
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